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How chaos is useful in the brain information processing is greatly unknown. Here, 
we show that the statistical property of chaos such as invariant measures naturally 
organized under a great number of iterations of chaotic mappings can be used for some 
complex computations, while the precise information of initial conditions which van- 
ishes in the course of iterations deos not matter for this kind of computations. The key 
observation of the present study is that computation using ergordicity of dynamical 
systems can be thought of as massively parallel Monte Carlo simulations. Here, to 
avoid difficulty in elucidating the ergordicity of dynamical systems, we propose com- 
putational schemes using the generalized class of one-dimensional chaos with explicit 
invariant measures. The validity of our results which connect chaos with parallel com- 
putation is checked by the precision computations of some transcendental numbers like 

7T. 
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One of the most important theoretical challenges in our brain information 
processing is how a great number of noisy neurons realize the parallelism of com- 
putation in the brain. Thus, in relation to this problem, it is of great interest to 
consider whether or how a kind of chaotic dynamics can be harnessed in brain 
information processing mechanisms. However, most studies concerning chaos in 
the brain are based on the empirical computer simulations on some specific con- 
ditions. The purpose of the present paper is to show that we can have some ideal 
model of parallel computation using chaos by purely analytical treatment. 

Let us consider the following deterministic process x n+ i = F(x n ), where 
x n G M . We can say that this dynamical system has the ergordicity if the 
following relation holds: 

1 N f 

j 1 ^ ^ Q( x i) = J M Q(x)v(dx), (1) 

2 = 1 

where p(dx) is the invariant measures on M and Q(x) is a regular function 
which makes R.H.S. of Eq. ([I]) finite. Essentially, Monte Carlo simulations can 
be regarded as evaluating R.H.S. of Eq. ([!]) by approximately calculating L.H.S. 
of Eq. ([!]) under the assumption of the ergordicity. However, it is difficult to get 
an explicit invariant measure fi(dx) for a given dynamical system and in general, 
we cannot tell whether the ergordicity holds or not for the system. Without 
explicit invariant measures, it is meaningless to consider the Monte Carlo type 
computations because the integral as R.H.S. of Eq. (|1|) is unclear. 

Quite recently, we discovered the two-parameters class of chaotic dynamical 
systems with explicit invariant measures^, ||. We call this class exactly solvable 
chaos because we can obtain the invariant measures as well as the general solu- 
tions analytically. Here, let us briefly review the models of exactly solvable chaos 
on the unit interval M = [0, 1] = I. The first example of exactly solvable chaos 
was given by Ulam and Von Neumann in 1947 as follows: 

= F(x n ) = 4x n (l x n ), (2) 

where the invariant measure has an explicit expression 

fi(dx) = p(x)dx = — (3) 
nJx(l — x) 

Next, in 1985, Katsura and Fukuda|]TJ generalized the Ulam=Neumann map to 
the following model: 

for < k < 1. Its invariant measure can be calculated as 

u(dx) = p(x)dx = , =, (5) 

2K{k)^x{l-x){l-k 2 x) 
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where K(k) is given by the complete elliptic integral 

K(k) =\ \ dU (6) 
^(l-u 2 )(l-k 2 u 2 ) 

The Katsura=Fukuda map for k = corresponds to the Ulam=Neumann map. 
Furthermore, the present author generalized the Katsura=Fukuda to the follow- 
ing model M: 

4x(l — x)(l — lx)(l — mx) 

i ( x) = 

1 — 2(1 + m + lm)x 2 + 8lmx 3 + (I 2 + m 2 — 2lm — 2l 2 m — 2lm 2 + l 2 m 2 )x^ 

(7) 

with two parameters m and / such that — oo < m < I < 1. If we set m = 0, 
we obtain the Katsura=Fukuda map. The explicit invariant measure of the 
generalized map has the form|4]] of 

fi(dx) = p(x)dx = . , (8) 

2K(l,m)Jx(l — x)(l — lx)(l — mx) 

where K(l,m) is given by the integral 

K(l, m) = f 1 dU (9) 

(I - u 2 )(l - lu 2 )(l - mu 2 ) 

Here, we call (^) the generalized Ulam=Neumann map. Furthermore, it is shown 
that the generalized Ulam=Neumann map (0) is the most generalized model 
of exactly solvable chaos on I, whose Lyapunov exponent is equal to log 2, i.e., 
we cannot further generalize the map (|?|). Thus, we can say that the class of the 
generalized Ulam=Neumann maps is universal in representing exactly solvable 
chaos on I. We note here that their dynamical zeta functions of the generalized 
Ulam=Neumann maps is the same as the one of the original Ulam=Neumann 
map||], by which it means that the statistical characters such as Kolmogorov= Sinai 
entropy, Lyapunov exponent, and topological entropy do not depend on the pa- 
rameters / and m. Figure 1 shows the shapes of the generalized Ulam=Neumann 
maps with different parameters m for I = 0.5. Interestingly, their shapes for m = 
—20, • • • , —100 are very similar to the Poincare plots of Belousov=Zhabotinski 
chemical reaction. 

Using the explicit expression of the invariant measures p(x) of the generalized 
Ulam=Neumann maps in Eq. (|8]), we obtain the following formula: 

1 N ri 1 

& n g G(x n ) = j o G(x)p(x)dx = (10) 



where G(x) = yx(l — x)(l — lx)(l — mx). K(l, m) is known to be a transcen- 
dental real number like n. Let us consider the problem of parallel computing 
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of K(l,m) using this model of exactly solvable chaos. We consider P indepen- 
dent processors labeled by j(l < j < P). First, we give an initial condition 
x\j e I to each processor. The initial conditions are different each other. Then 
each processor iterates the map (|7|) as x n+ xj = F(x n j) to compute the number 
= jjJ2iLiG(x n) j). Finally, after averaging the above numbers, we have 

an estimated value of K(l,m) as 

2K(l,m) ^1^2%,-) « 2K(l,m). (11) 

To confirm the validity of this proposed parallel algorithm, we computed 2K(l, m) 
for (l,m) = (0,0), (0.5, -100) and (0.5,-200). We set that P = 100, N = 10 6 
and xij = 0.0073333371 j for j = 1, • • • , 100. 



Table 1: Results for P = 100 and N = 1000000 





(/,m) = (0,0) 


(l,m) = (0.5,-100) 


(l,m) = (0.5,-200) 


2K{l,m) 


3.14167445904031 


0.804839363844501 


0.618717074365706 


2K{l,m) 


3.14159265358979 


0.804826090564245 


0.618817767268870 



Table 1 shows that our numerical estimated values 2K(l,m) using the stan- 
dard double precision representations of real numbers neared the exact values 
2K(l, m) as \2K{1, m) - 2K(l, m) \ < 0.0002. 

Furthermore, as is shown in Fig.2, our estimated values 2K(l,m) for (l,m) = 
(0.5, 0.3) using exactly solvable chaos converge faster than those using the usual 
Monte Carlo method to evaluate integrals on I, although both of the mean square 
deviations are inversely proportional to the iteration number N. This fastness 
of convergence of our method over the Monte Carlo method is owing to the fact 
that G(x) = \/x(l — x)(l — lx)(l — mx) < 1 < , 1 = for x,y E I 

and < m < I < 1. 

This work was supported by the special researcher's program to promote ba- 
sic sciences at RIKEN. 
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Fig. 1: The generalized Ulam=Neumann maps (7) for (l,m) 
(0.5,0), (0.5, -10), •••,(0.5, -100). 
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Fig. 2: Mean squared deviations defined as -p J2f=i{^^{ x i,j) ~ 2_ftT(Z, m)} 2 with 
(I, m) = (0.5, 0.3) versus the number of iterations N for the present method 
(Chaos Computing) and for the usual method (Monte Carlo), where the uni- 
form random numbers were generated by the function ran(seed=564789 11) in 
the FORTRAN77 library to determine the initial conditions x±j for the Chaos 
Computing and the Monte Carlo paths. 
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Generalized Ulam=Neumann Maps for 1=0.5 




Chaos Computing versus Monte Carlo for (l,m,P)=(0.5,0.3,100) 
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